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Plasmon and polariton modes are derived for an ideal semi-infinite (half-space) 
plasma and an ideal plasma slab by using a general, unifying procedure, based on 
equations of motion. Maxwell's equations and suitable boundary conditions. Known 
results are re-obtained in much a more direct manner and new ones are derived. The 
approach consists of representing the charge disturbances by a displacement field in 
the positions of the moving particles (electrons). The dielectric response and the elec- 
tron energy loss are computed. The surface contribution to the energy loss exhibits 
an oscillatory behaviour in the transient regime near the surfaces. The propagation 
of an electromagnetic wave in these plasmas is treated by using the retarded elec- 
tromagnetic potentials. The resulting integral equations are solved and the reflected 
and refracted waves are computed, as well as the reflection coefficient. For the slab 
we compute also the transmitted wave and the transmission coefficient. Generalized 
Fresnel's relations are thereby obtained for any incidence angle and polarization. Bulk 
and surface plasmon-polariton modes are identified. As it is well known, the field in- 
side the plasma is either damped (evanescent) or propagating (transparency regime), 
and the reflection coefficient for a semi-infinite plasma exhibits an abrupt enhance- 
ment on passing from the propagating regime to the damped one (total reflection). 
Similarly, apart from characteristic oscillations, the reflection and transmission coef- 
ficients for a plasma slab exhibit an appreciable enhancement in the damped regime. 
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I. INTRODUCTION 

After the discovery of bulk plasmons in an infinite electron plasma, there was a great 
deal of interest in plasmons occurring in structures with special geometries, like a half-space 
(semi-infinite) plasma, a plasma slab of finite thickness, a two-plasmas interface (two plas- 
mas bounding each other), a two-dimesional sheet with an aperture, a slab with a cilindrical 
hole, structures with surface gratings or regular holes patterns, layered films, cilindrical 
rods and spherical particles, etc. There is a vast literature on various structures with 
special geometries exhibiting plasmon modes. These studies were aimed mainly at iden- 



tifying new plasmon modes, like the surface plasmons, J|- JjJ accounting for the electron 
energy loss experiments and exploring the interaction of the electron plasma with electro- 

nn 

magnetic radiation (polariton excitations). |12l|-|24l| More recently, a possible enhancement of 
the electromagnetic radiation scattered on electron plasmas with special geometries enjoyed 



a particular interest. 
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ji all these studies the plasmon and polariton modes are of 

n n 

fundamental importance. [28!|-[32| The methods used in deriving such results are of great 
diversity, resorting often to particular assumptions, such that the basic underlying mecha- 
nism of plasmons or polaritons' occurrence is often obscured. The need is therefore felt of 
having a general, unifying procedure for deriving plasmon and polariton modes in structures 
with special geometries, as based on the equation of motion of the charge density. Maxwell's 
equations and the corresponding boundary conditions. Such a procedure is presented in this 
paper for an ideal semi-infinite plasma and an ideal plasma slab. 

We represent the charge disturbances as 6n = —ndivu, where n is the (constant, uniform) 
charge concentration and u is a displacement field of the mobile charges (electrons). This 
representation is valid for Ku(K) -C 1, where K is the wavevector and u(K) is the Fourier 
component of the displacement field. We assume a rigid neutralizing background of positive 
charge, as in the well-known jellium model. In the static limit, i.e. for Coulomb interaction, 
the lagrangian of the electrons can be written as 

L = J dr ^mnu^-^J dr'U{\r - r'\)6n{r)6n{r') +ej dr^{r)6n{r) , (1) 



where m is the electron mass, U{r) = e^/r is the Coulomb energy, — e is the electron charge 
and $(r) is an external scalar potential. Equation ([T]) leads to the equation of motion 

mil = ngrad J dr'U{\r — r'\)divu(r') + egrad^, (2) 

which is the starting equation of our approach. We leave aside the dissipation effects (which 
can easily be included in equation ([2])). 

By using the Fourier transform for an infinite plasma it is easy to see that the eigenmode of 
the homogeneous equation ^ is the well-known bulk plasmon mode given by Up = Anne^ /m. 
On the other side, equation 5n = —ndivu is equivalent with Maxwell's equation divEi = 
—4:7ie6n, where Ej = Anneu is the internal electric field (equal to — 47rP, where P is the 
polarization). Making use of the electric displacement D = —grad^ = £(D + Ej), we get the 
well-known dielectric function e = 1 — ujp/u'^ in the long- wavelength limit from the solution 
of the inhomogeneous equation ([2]). Similarly, since the current density is j = —enu, we get 
the well-known electrical conductivity a = iUp/Anuj. 

We apply this approach to a semi-infinite plasma and a plasma slab. First, we derive the 
surface and bulk plasmon modes and obtain the dielectric response and the electron en- 
ergy loss for a semi-infinite plasma. The surface contribution to the energy loss exhibits 
an oscillatory behaviour in the transient regime near the surface. Further on, we consider 
the interaction of the semi-infinite plasma with the electromagnetic field, as described by 
the usual term (1/c) / c/rjA — / drp^ in the lagrangian, where A is the vector potential, 
p = endivu is the charge density and $ is the scalar potential. We limit ourselves to the 
interaction with the electric field, and compute the refiected and refracted waves, as well 
as the reflection coefficient. Generalized Fresnel's relations are obtained for any incidence 
angle and polarization. We find it more convenient to use the radiation formulae for the 
retarded potentials, instead of using directly the Maxwell's equations, and the resulting in- 
tegral equations are solved. Bulk and surface plasmon-polariton modes are identified. The 
field inside the plasma is either damped (evanescent) or propagating (transparency regime), 
and the reflection coefficient exhibits an abrupt enhancement on passing from the propagat- 
ing to the damping regime (total refiection). Finally, we give similar results for a plasma 
slab, where we compute also the transmitted field and the transmission coefficient. Apart 
from characteristic oscillations, the refiection and transmission coefficients for a plasma slab 



exhibit an appreciable enhancement in the damped regime. The present approach can be 
extended to various other plasma structures with special geometries. 



II. PLASMA EIGENMODES 



We consider an ideal semi-infinite plasma extending over the half-space z > (and bounded 
by the vacuum for z < 0). The displacement field u is then represented as (v, U3)9{z), where 
V is the displacement component in the (x, ?/)-plane, M3 is the displacement component along 
the z-direction and 9{z) = 1 for z > and 9{z) = for z < is the step function. In equation 
of motion ((21) divu is then replaced by 

divu = (^divY + e{z) + U3{0)6{z) , (3) 
where ^3(0) = u^i^r, z = 0), r being the in-plane (x, y) position vector. Equation (l2|) becomes 



mil = ne^grad J dr'dz' —, = divvir' .z') + 

Y (r— r')2+(z— 2;')2 L 



+ne'^grad J dr'—. — = — ='U3(r', 0) + egrad^ 

W (r— r')2+22 



+ 

(4) 



for 2; > 0. One can see the (de)-polarizing field occurring at the free surface z = (the 
second integral in equation (111)). 

We use Fourier transforms of the type 

u(r, z;t) = Y, f duju{k, z; cu)e*^'"e-*^* (5) 
k 

(for in-plane unit area), as well as the Fourier representation 

1 9-7r 



\/r^ + 2;2 ^ k 

for the Coulomb potential. Then, it is easy to see that equation ^ leads to the integral 
equation 



uj'^v = - 
2 



1 2 f°° 1 I -k\z-z'\ 1 2 f°° 7 d _k\z-z'\ ^e/c^ . . 



and iku^ = where we have dropped out for simplicity the arguments k, z and cu. The v- 
component of the displacement field is directed along the wavevector k (in-plane longitudinal 



waves). This integral equation can easily be solved. Integrating by parts in its rhs we get 



uj\ = u^v - lulvoe - — $ , (8) 
^ 2 ^ m 



hence 

(9) 

^ f^kuip $g e~^-^ e 

where vq = v{z = 0)^ = ^{z = 0) and $' = One can see the surface contributions 
(terms proportional to $0^^''^) and bulk contributions ($, $'-terms). 

The solutions given by equations ([9]) exhibit two eigenmodes, the bulk plasmon = ojp and 
the surface plasmon ujs = ujp/^/2, as it is well known. Indeed, the homogeneous equation 
(l8|) ($ = 0) has two solutions: the surface plasmon v = VqC'^^ for u"^ = ujp/2 and the bulk 
plasmon Vq = for u'^ = ^p- Making use of this observation we can represent the general 
solution as an eigenmodes series 



for z > 0, where T;(k, —k) = — f (k, k), and ikusCk, z) = Then, it is easy to see that 

the hamiltonian H = T + U corresponding to the lagrangian L = T — U given by equation 
([I]) becomes 

T = nni X^k i'oO^)i'oO^) + n-m X^k^ K)v(k, k) 

(11) 

U = 2Tm^e^ Ek ^o(k)t;o(k) + Ann^e^ Ek« ^^*(k, K)f (k, k) , 
where T is the kinetic energy and U is the potential energy. We can see that this hamiltonian 
corresponds to harmonic oscillators with frequencies tOs = ^^p/ V2 and tOb = 



Making use of Ej = 47meu and equations ([9]) we can write down the internal field (polar- 
ization) as 

-C/Xl^', - (^2_^2)(2<^2_^2)e ^2_^2 



(12) 



IZ/yl^K, - (a;2_^2)(2a;2_<^2)e ^2_^2 



where E± is directed along the in-plane wavevector k and E\\ is parallel with the 2;-axis 
(perpendicular to the surface z = 0). This is the dielectric response of the semi-infinite 
plasma to an external potential. 



We take an external potential of the form $(k, 2;) = <l>°(k)e*'^^ (leaving aside the frequency 
argument w), and get the electric displacement D±(k,z) = — zk$°(k)e*''^ and D^\(k, z) = 
— m$°(k)e*''^ from D = —grad^. We can see that the surface terms do not contribute to 
this response, as expected, since these terms are localized. Making use of Ej = (1/e — 1)D, 
we get the well-known dielectric function e{K,Lj) = 1 — oj^juJ^ in the long-wavelength limit. 



III. ELECTRON ENERGY LOSS 



It is well known that the energy loss per unit time (stopping power) is given by 



for an electron moving with velocity v = (vx,t>||), where the field E, is taken at r = v^t 
and z = for t > {z > 0). It is assumed that the electron energy is sufficiently 
large and the energy loss is small enough to use a constant v in estimating the rhs of 
equation (fT3l) . The potential created by the electron is given by the Poisson equation 
A$ = 47re(5(r — \j_t)6{z — v\\t), whence, by making use of the Fourier representation (I6l), we 
get 

•^(k, ^; 00) = -, ^^f-TT^e-^C^---)^/''" . (14) 

[UJ — kV_|_)^ + K^f y 

We introduce this potential in equations (fT2ll and compute the energy loss given by equation 
( fT3l) . It contains two contributions, one associated with the bulk plasmons, 

iuj 2v\\ 

and another arising from surface effects, 



2 2 

e 



eV y / du ^ ^;||(.kv^-^|) (w,-.), . 

In performing the w-integrations in equations (fTSi) and (fT6l) we retain only the plasmon 
contributions arising from the poles uj = ojp and uj = iJp/^/2. For normal incidence {v± = 0, 
f|| = v) we get easily the well-known bulk contribution Pf, = (^—e'^oOp/vj \n{vko/ujp), where 
ko is an upper cut-off (associated, as usually, with the ionization energy, or with the inverse 
of the mean inter-particle spacing, etc), and the surface contribution 

Ps = ^ (^■\/2 sin u pt / V2 — sinupt^ . (17) 




Figure 1: Electromagnetic plane wave Eg, with wavevector K, incident on the surface z = 0. 

We can see in equation (fT71) the oscillatory behaviour of the stopping povi^er arising from the 
surface effects in the transient regime near the surface. 

IV. INTERACTION WITH THE ELECTROMAGNETIC FIELD. POLARITONS 

We assume a plane v^^ave incident on the plasma surface under angle a. Its frequency is 
given hj u = cK, where c is the velocity of light and the wavevector K = (k, k) has the 
in-plane component k and the perpendicular-to-plane component k, such as k = i^sina 
and K = Kcosa. In addition, k = k{cosip,smip). The electric field is taken as Eq = 
i?o(cos P, 0, — sin /5)e*'''"e*'*^e~*'^*, and we impose the condition cos P sin a cos — sin (3 cos a = 
(transversality condition KEq = 0). The angle /? defines the direction of the polarization 
of the incident field. The geometry of the incident wave is shown in Fig. 1. 

In the presence of an electromagnetic wave we use the equation of motion 

u^u = — E + —Eoe'"' , (18) 
m m 

for z > 0, where E is the polarizing field; in equation (fTSl) we have preseved explicitly only 
the 2;-dependence {i.e. we leave aside the factors e*'''"e~*'^*). We find it convenient to employ 



the vector potential 



A(r, z;t) = lJ dr' J d.- ^^^ ' ^ ^^^^ (19) 



and the scalar potential 

t(r..;i) = /dr-/d.' ^''-'--;-^^-> . (20) 

where j = —neu6{z)e^^^e~^'^^ is the current density, p = nedivu = 

ne (ikv + 9{z)e^^''e~^'^*^ + neu3{0)6{z)e^^''e~'^'^^ is the charge density and 



R = {r — r')^+ {z — z'Y . The integrals in equations (fTOl) and (l20l) implies the 
known integral(33| 

/ rfx Jo f A; Vx2 - e'""/'^ = -e*"l^l , (21) 

where Jq is the zeroth-order Bessel function of the first kind (and tu^/c^ = + fc^). It 
is convenient to use the projections of the in-plane displacement field v on the vectors k 
and = A;(— sin cosy)), k^k = 0. We denote these components by v\ = kv/k and 
t>2 = k±\/k, and use also the components Ei = kE/A;, E2 = k^Fi/k and similar ones for 
the external field Eq. We give here the components of the external field 

Eqi = Eq cos P COS ip , £'02 = —Eq COS P sinip , Eo3 = —Eq sin (3 . (22) 

One can check immediately the transversality condition E^ik + Eq^.^ = 0. Making use of 
E = — — grad^, equations (fTOl) and (l20l ) give the electric field 

El = -27iineK jQdz'vi{z')e'^\'-^'\ - 2Tme^ jQdz'u^{z')^e'''\'-''\ 



E2 = -2TTine^ /o dz'v2iz')e'''^'-''^ (23) 

E3 = 27me^ ^Qdz'vi{z')-^e^'^\'^^'''\ - 2-Kine^ jQdz'u3{z')e^'^^'''''''^ + Aimeu^ 

for z > 0. It is worth observing in deriving these equations the non-intervertibility of the 
derivatives and the integrals, according to the identity 

— dzf{z)—e dzf{z)e -2tKf{z) (24) 







dz Jo dz 

for any function /(2;), z > 0; it is due to the discontinuity in the derivative of the function 
for z = z . Now, we employ equation of motion (fTSl) in equations (l23l) and get the 



IK 

e 



integral equations 



uj'v2 = /o dz'v2iz')e'''\^-^'\ + ^Eo2e'-^ (25) 



^'^3 = # /o dz'v,{z')£e''^\^-^'^ - ^ /o dz'us{z')e'^\^-^'\ + u;>3 + ^Eosc'^' 
for the coordinates fi 2 and M3 in the region z > 0. 

The second equation (l25ll can be solved straightforwardly by noticing that 

^ / dz'v2iz')e'''\'-''\ = -K^ I dz'v2{z')e'''\'-''\ +21KV2 . (26) 
oz^ Jo Jo 

We get 



d^V2 



Q^2 + ('^ - ^p/c = . (27) 
The solution of this equation is 



Op Enn K[K — K 

V2 = '-^- \, ' , (28) 

where 

K = \J n'^ — i^p/c^ = --y/cu^ cos^ a — uj'^ . (29) 

The wavevector /t'can also be written in a more familiar form k = {u / c)\/ e — sin^ a, where 
£ = 1 — cUp/u;^ is the dielectric function. The corresponding component of the (total) electric 
field (the refracted field), can be obtained from equation (fT8l ): it is given by (mtu^/e) V2. For 
< ojpjc^ (u;cosa < Up) this field does not propagate. For > oj^/c^ {uj greater than the 
transparency edge Up/ cos a) it represents a refracted wave (transparency regime) with the 
refraction angle a given by Snell's law 

sin a' 1 ^ I r- ,„ ^ 
= , = • 30 

The polariton frequency is given by 

C^2 = c^if 2 = ^2 + c^K'^ ^ (31) 

as it is well known, where K"^ = k"^ + k"^. 
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a) a=7c/6 

b) a=7i/5 

c) a=7i/4 

d) a=Ji/3 

e) a=7t/2 




Figure 2: Reflection coefficient for a semi-infinite plasma for /3 = vr/G and various incidence angles 
a. One can see the shoulder occurring at the transparency edge ujp/ cos a and the zero occurring 
at = t^2^ (1 _ tan^ q) for a = /5 = n/Q {R2 = 0,(/7 = 0). 



The first and the third equations (12511 can be solved by using an equation similar vv^ith 
equation ( !26l ) and by noticing that they imply 

dvi 



K U3 = ik- 



dz 



We get 



and 



2eE, 



Vl 



Us 



01 



2eE, 



03 



K,' (k, — K,'] f 

K\K — K] / 

V / „IK z 



(32) 



(33) 



(34) 



Similarly, the corresponding components of the refracted field are given by equation (fTSll . It 
is easy to check the transversality condition V\k ^ u^k' = (and the vanishing of the bulk 
charge ne (^zkv + = 0). 

We can see that the polarization field E in equation (fT8l ) cancels out the original incident 



field Eo and gives the total, refracted field mcu^u/e inside the plasma. This is an illustration 
of the so-called Ewald-Oseen extinction theorem. (itI. Is^ 



It is worth investigating the eigenvalues of the homogeneous system of integral equations 



(l25ll . for parameter k given by k = yuj'^/c^ — k'^. Such eigenvalues are given by the roots of 
the vanishing denominator in equations (l33l ) and (l34l) . i.e. by equation kk + k"^ = 0. This 
equation has real roots for u only for the damped regime, i.e. for k = i |k| and k = i k . 
Providing these conditions are satisfied, there is only one acceptable branch of excitations, 
given by 

= " , . 35 



We can see that a; ~ cfc in the long wavelength limit and it approaches the surface-plasmon 

frequency u ~ ujp/\/2 in the non-retarded limit {ck oo). These excitations are surface 

/ 

— K Z 

plasmon-polariton modes. We note that they imply f 2 = and Wi, % ~ e .In addition, 
a careful analysis of the homogeneous system of equations (l25l) reveals another branch of 
excitations, given hy ui = oOp, which, occurring in this context, may be termed the bulk 
plasmon-polariton modes. They are characterized by f 2 = and fi(k, 0) = 0. For all these 
modes we have M3 = ic^k/ (uj'^ — c^k"^ — uj^ 



dz ■ 



In order to get the reflected wave (the region 2; < 0) we turn to equations (i23ll and use therein 
the solutions given above for f 1 2 and M3. It is worth noting here that the discontinuity term 
ciJpUs does not appear anymore in these equations (because z > and z < and we cannot 
have z = z'). The integrations in equations (l23ll are straightforward and we get the held 

E, = Eoi^^ ■ , ,, e--"^ , 36 

K + K KK + k"^ 

E2 = Eo2^^e-'^^ (37) 

and 

' ' ; 2 

K + K KK + k-^ 

We can see that this field represents the reflected wave {k — > —k), and we can check its 
transversality to the propagation wavevector. Making use of the reflected field E^e// given 
by equations fl36l) -fl38l) and the refracted field 'Ej-efr obtained from equations (fTSll and (l23l) 
(Erefr = E + Eq = mcu^u/c) ouc cau check the continuity of the electric field and electric 
displacement at the surface {z = 0) in the form Ei^2refi + -^01,2 = -E'i,2re/r, E^refi + -E'oa = 
eEsrefr, whcrc e = 1 — Up/uj'^. The angle of total polarization (Brewster's angle) is given 
by KK — k'^ = 0, or tan^ a = 1 — tOp/to"^ = £ (for a < tc/A). The above equations provide 



generalized Fresnel's relations between the amplitudes of the reflected, refracted and incident 
waves at the surface for any incidence angle and polarization. They can also be written by 
using u;^ = tUp/ (1 — e), where e is the dielectric function. 

The reflection coefficient R = {Eirefif / {Eof can be obtained straightforwardly from the 
reflected fields given by equations (i36ll-(l38ll. It can be written as 



R = Ri cos p sin (p + R2 (cos (3 cos ^9 + sin 



(39) 



where 



and 



Ri 



uj'^ cos^ a 



LUp — uj cos a 



uj^ cos^ a 



ujp + uj cos a 



R2 



cos awtt;^ cos^ a 



Up — u sin^ a 



cosa^/uj'^ cos^ a — ujp + u sin^ a 



(40) 



(41) 



The first term in the rhs of equation (l39l) corresponds to P = {ip = n/2; s-wave, electric 
field perpendicular to the plane of incidence, while the second term corresponds to P = a 
[if = 0; p-wave, electric field in the plane of incidence). It is easy to see that there exists a 
cusp (shoulder) in the behaviour of the function Rioj), occurring at the transparency edge 
ui = Up/ cos a, where the refiection coefficient exhibits a sudden enhancement on passing from 
the propagating regime to the damped one, as expected (total refiection). The condition 
for total refiection can also be written as sin a = ^/e, where R = 1 (-Ri,2 = 1), as it is 
well known. For illustration, the refiection coefficient is shown in Fig. 2 for /? = 7r/6 and 
various incidence angles. The refiection coefficient is vanishing at u"^ = uj^j (1 — tan^ a) for 
a = /? < 7r/4 (i?2 = 0, = 0). 



V. PLASMA SLAB 



We consider an ideal plasma slab of thickness (i, extending over the region < 2; < 
d and bounded by the vacuum. The displacement field u can be represented as 
(y ^uz)\§{z) — Q{z — d)\, where v is the displacement component in the (x, |/)-plane and 
U3 is the displacement component along the z-direction. The approach presented above for 
a semi-infinite plasma can easily be extended to this case. The analogous of the equation of 
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Figure 3: Reflection coefficient for a slab of thickness d {dtOp/c = 1) for /3 = 0, = 7r/2 (s-wave) 
and a few incidence angles a. Its slope is continuous at the transparency edge (wcosa = cOp). The 
oscillations occurring in the transparency regime are too small to be visible in Figure. 

motion (111) exhibits now two polarization contributions, arising from the two surfaces. The 
dielectric response similar to equation ([9]) is given by 

(42) 

and ikus = |f , where = ^{z = 0), $rf = $(z = d), < z < d. The electric field is given 
by E± = innev and E\\ = Ajmeu^. One can see that, beside the bulk plasmon mode tu^, 
there appears two surface modes given by ujI (l ± e'^'^) /2, as it is well known. For d ^ oo 
equation ( l42ll becomes the first equation ([9]) for the semi-infinite plasma. For d — > we get 
the well-known plasma frequency ^J {T.nnse^ /m) k for a sheet with surface electron density 
Us = nd. 

The bulk contribution to the energy loss is the same as for the semi-infinite plasma. We 
compute the surface contrbution to the electron energy loss for kd ^ Upd/v ^ 1, i.e. for a 
fast electron moving with velocity v, which, however, spends enough time in the sample to 
excite plasmons. For normal incidence the surface contribution consists of two oscillatory 
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Figure 4: Reflection coefficient for a slab of thickness d {dojp/c = 1) for a = /3, 93 = (p-wave) and 
a few incidence angles a. It exhibits a local maximum {R = 1) for uj = ujp and small oscillations in 
the transparency region uj cos a > tOp (too small to be visible in Figure). In addition, it is vanishing 
for (jj"^ = ujp/ (1 — tan^ a), a < 7r/4, as one can see in Figure for a = 7r/6 (curve a). 

terms 

Ps = —^-^ (^V2smujpt/y/2 — sinujptj — 

(43) 

\^smujp (d/v — t) /^/2 — sintOp (d/v — t) , 
corresponding to the two surfaces, for < t < d/v. The total energy loss during the passage 
through the slab is given by 

dtPs - dtPs = -TT ( V2 - 1) ^ . (44) 

We use again the equation of motion (fTSll and the retarded potentials given by equations 
(fT9|) and (l20|) in order to get the refracted field (field inside the slab), refiected {z < 0) and 
transmitted {z > d) fields. The polarization field is given by the same equations ( l23l ). where 
the z-integration is limited to the region Q < z < d. The same holds for the equations of 
motion (i25l l. We solve these equations by the same method used above. Within the slab 
we have two waves of the form e^*** ^, one being the refracted wave through the first surface 
(z = 0), the other being the reflected wave on the second surface {z = d). The wavevector k 




a) oc=7i/6 

b) a=7t/4 

c) a=Ji/3 



d—vt 



is given by the same equation (l29l) . and the transparency edge is given by the same condition 
u) cos a = as for a semi-infinite plasma. We get 



V2 = A-i 



gift z _____^2iK a ^ 



K + K, 



where 



and 



where 



A, 



2eEo2 









vi = Ai 



K + K KK + k'^ 



Ai 



2eE, 



01 



K \ K — K ] \ K + K 



/\2 



KK + k'^ 



(45) 



(46) 



(47) 



(48) 



mw^ {k + kY {kk + kY - (/« - n'y {kk - k^Y e2« 
the third component can be obtained from k'^us = ik (dvi/dz). One can check the transver- 
sality of these waves and can compute the dispersion relations for the eigenvalues (bulk and 
surface plasmon-polaritons) in the like manner as for the semi-infinite plasma. 



The refiected field is given by 

El = Eoi (l - e^-' 



(k+k') (kk'+Zc^) — (k— k') (kk'— fc^) e^*" 



(49) 



2iK d 
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E2 = Eq2 1^1 — e 
and E3 = -Eo3 {Ell Em). 

>From the above results one can check the continuity of the electric field and electric dis- 
placement as well as the angle of total polarization given by tan^ a = \ — uj^joj^ = e. If we 
take formally e^*'' — > we recover all the fields for the semi-infinite plasma. Indeed, for the 
semi-infinite plasma all the integrations to 2; — cx) are taken by assuming a vanishing factor 



-/IZ 



fi > 0, and letting /i go to zero. If we preserve this factor for the slab, it gives rise 
to factors of the form e^*** '^e~'^'^, which are vanishing for d 00. The limit d — > (plasma 
sheet) cannot be taken directly on the above results {ujp ~ 1/v^, k! ~ i^^p/c), because of 
the discontinuities arising from the 6'-function. The calculations for a plasma sheet with a 
finite (superficial) charge density must be done separately. They are left, together with 
other related results, for a forthcoming publication. The limit k d -C 1 {Kd -C 1) can be 
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Figure 5: Transmission coefficient for a slab of thickness d {dtOp/c = 1) for /3 = 0, (/9 = 7r/2 (s- 
wave) and a few incidence angles a. One can see the characteristic cusp at the transparency edge 
Lj COS a = iVp and the peak ocurring below this edge. The oscillations occurring in the transparency 
regime are too small to be visible in Figure. 



taken directly on the formulae given here. It corresponds to wavelengths much longer than 
the thickness of the slab. 

The reflection coefficient for the plasma slab R = |Ere/i|^/ |Eo|^, where the reflected field 
is given by equations (llQl) . has a different structure than the reflection coefficient for the 
semi-infinite plasma. It can be written as 

A , 2 



1-e 



2iK 



R\ cos^ /3 sin^ Lp + R2 {cos^ (3 cos^ (p + sin^ 



where 



and 
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K + K 



Ri 



(50) 



(51) 



(52) 



[k + k) [kk + fc2) — (k — k) [kk — k"^) 
The reflection coefficient given by equation (l50l) is shown in Figs. 3-4 for /? = 0, = 7r/2 
(s-wave) and, respectively, a = /3, = (p-wave) and dojp/c = 1. The reflection coefficient 
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Figure 6: Transmission coefficient for a slab of thickness d [dwp/c = 1) for a few incidence angles 
a = (3 and (/? = (p-wave). One can see the two peaks occurring below the transparency edge 
u! COS a = ujp (the cusp in Figure) and the zero for uj - 
transparency regime are too small to be visible in Figure. 



uOp. The oscillations occurring in the 



exhibits characteristic oscillations arising from the exponential factor in equations (!50l)- (i52ll 
and has an abrupt enhancement in the damping regime. In addition, i?2 is vanishing for 

o;^ = Up/ (1 — tan^ a) (a < 7r/4) and -R2 = 1 for = ujp. 

The transmitted field (region z > d) is given by 
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and E'a = Eq^ (Ei/Eqi). One can check the continuity of the electric field and electric 
displacement at the surface z = d. In the limit d ^ 00 the transmitted field is vanishing. 
The transmission coefficient given by T = |Et,.|^ / |Eo|^, where Ej^ is given by equations 



53ll. can be written as 



T = 16k'^ 



K 



Ri cos^ P sin^ if + 
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(54) 



where R12 are given by equations (ISTll and (l52l) . This transmission coefficient is shown 
in Figs. 5-6 for /? = 0, = 7r/2 (s-wave) and, respectively, a = P, ip = (p-wave) and 
dujp/c = 1. Beside the characteristic cusp occurring at the transparency edge (cjcosa = 
Up), the transmission coefficient exhibits an appreciable enhancement below this edge. For 
a = P, ip = (p-wave) and u = Up the reflection coefficient attains the value unity and 
the transmission coefficient vanishes. The fields derived above can be viewed as generalized 
Fresnel's relations for a plasma slab. 



VI. CONCLUSIONS 



The approach presented here is a quasi-classical one, valid for wavelengths much longer 
than the amplitude of the Fourier components of the displacement field u. This is not 
a particularly restrictive condition for the classical dynamics of the electromagnetic field 
interacting with matter. When this condition is violated, as, for instance, for wavelengths 
much shorter than the mean separation distance between electrons, there appear both higher- 
order terms in the equations of motion and the coupling to the individual motion of the 
electrons. These couplings affect in general the dispersion relations and introduce a finite 
lifetime (damping) for the plasmon and polariton modes. 

Making use of the equations of motion for the displacement field u and the radiation formu- 
lae for the electromagnetic potentials, we have computed herein the plasmon and polariton 
modes for an ideal semi-infinite electron plasma and an ideal plasma slab of finite thick- 
ness, as well as the dielectric response, the electron energy loss, the refiected and refracted 
waves and the refiection coefficient. For the semi-infinite plasma we have identified the bulk 
and surface plasmon-polariton modes and for the plasma slab we have computed also the 
transmitted wave and the transmission coefficient. It was shown that the stopping power 
due to the surface effects has a characteristic oscillatory behaviour in the transient regime 
near the surfaces. The field inside the plasma is either damped (evanescent) or propagat- 
ing, as it is well known, and the refiection coefficient for the semi-infinite plasma exhibits a 



sudden enhancement on passing from the propagating to the damped regime, as expected. 
The transparency edge is given by a; cos a — cOp, where a is the incidence angle, cu is the 
frequency of the incident wave and Up is the plasma frequency. Apart from characteris- 
tic oscillations, the reflection and transmission coefficients for the plasma slab exhibit an 
appreciable enhancement below the transparency edge. 

Other effects related to the dynamics of a semi-infinite electron plasma, or, in general, 
various plasmas with rectangular geometries, can be computed similarly by using the method 
presented here. The method can also be applied to plasmas with other, more particular, 
geometries. The dissipation can be introduced (as for metals) and a model can be formulated 
for dielectrics, amenable to the method presented here. This will allow the treatment of more 
realistic cases as well as various interfaces, in particular plasmas (or metals) bounded by 
dielectrics. These investigations are left for forthcoming publications. 
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